DILATATION OF A ONE-DIMENSIONAL NONLINEAR CRACK 
IMPACTED BY A PERIODIC ELASTIC WAVE* 
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Abstract. The interactions between linear elastic waves and a nonlinear crack with finite com- 
pressibility are studied in the one-dimensional context. Numerical studies on a hyperbolic model of 
contact with sinusoidal forcing have shown that the mean values of the scattered elastic displace- 
ments are discontinuous across the crack. The mean dilatation of the crack also increases with the 
amplitude of the forcing levels. The aim of the present theoretical study is to analyse these nonlin- 
ear processes under a larger range of nonlinear jump conditions. For this purpose, the problem is 
reduced to a nonlinear differential equation. The dependence of the periodic solution on the forcing 
amplitude is quantified under sinusoidal forcing conditions. Bounds for the mean, maximum and 
minimum values of the solution are presented. Lastly, periodic forcing with a null mean value is 
addressed. In that case, a result about the mean dilatation of the crack is obtained. 

Key words. Elastic wave scattering, contact acoustic nonlinearity, nonlinear jump conditions, 
periodic solutions of differential equations, dependence of solutions on parameters 

AMS subject classifications. 34C11, 35B30, 37C60, 74J20 
1. Introduction. 

1.1. Aims. The modeling of interactions between ultrasonic waves and cracks 
is of great interest in many fields of applied mechanics. When the wavelengths are 
much larger than the width of the cracks, the latter are usually replaced by zero- 
thickness interfaces with appropriate jump conditions. Linear models for crack-face 
interactions have been widely developed [21j . However, these models do not prevent 
the non-physical penetration of crack faces. In addition, laboratory experiments have 
shown that ultrasonic methods based on linear models often fail to detect partially 
closed cracks |25j . 

A well-known nonlinear model for cracks is the unilateral contact model [22l[T2]. A 
more realistic hyperbolic model accounting for the finite compressibility of crack faces 
under normal loading conditions has been presented for applications in engineering [2] 
and geomechanical contexts [4. The well-posedness and the numerical modeling of 
the latter model in the context of 1-D linear elastodynamics was previously studied in 
[15] . The generation of harmonics and the distortion of simulated scattered velocity 
and stress waves were also observed. 

Subsequent numerical experiments have brought to light an interesting property. 
In the case of a sinusoidal incident wave, the mean values of the elastic displacements 
around a crack were found to be discontinuous. The simulations conducted also indi- 
cated that the mean dilatation of the crack increases with the amplitude of the forcing 
levels. This purely nonlinear phenomenon (the linear models predict no dilatation) is 
of physical interest: experimenters can measure the dilatation |14j and use the data 
obtained to deduce the properties of the crack. 

A preliminary theoretical study based on the use of a perturbation method was 
presented in 16J. An analytical expression for the dilatation was obtained in that 
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study and successfully compared with the numerical results. However, this expres- 
sion gives only a local estimate, valid in the particular case of the hyperbolic model 
subjected to very low sinusoidal forcing levels. The aim of the present study is to 
provide a theoretical analysis which can be applied to any range of forcing levels and 
to a larger number of contact models. 

1.2. Design of the study. The present paper is organized as follows: 

• In section 2, the physical configuration is described. The family of strictly 
increasing convex jump conditions dealt with in this study is introduced. 
Numerical simulations of the process of interest are presented; 

• In section 3, the problem is reduced to a nonautonomous differential equation. 
Up to section 6, only the case of sinusoidal forcing, which allows a complete 
understanding of the phenomena involved, is studied; 

• In section 4, preliminary results on the differential equation are presented. 
Classical tools for dynamical systems are used: Poincare map, phase portrait, 
lower and upper solutions [10l [18] ; 

• In section 5, the main qualitative results of this study are presented. The 
mean, maximum and minimum aperture of the crack are bounded, and local 
estimates for small forcing are also proved; 

• In section 6, some of the previous results are extended to non-monochromatic 
periodic forcing conditions. The increase of the mean dilatation with the 
forcing parameter is proved in theorem 16.31 The most general quantitative 
result of the paper is given in equation (|6.6p : 

• In section 7, some conclusions are drawn about the physical observables. 
Some future perspectives are also suggested. 

2. Statement of the problem. 




R 



Fig. 2.1. Elastic media Qq and f2i with rough contact surfaces, under constant static stress p. 
Static (left) and dynamic (right) case, with incident (I), reflected (R) and transmitted (T) waves. 



2.1. Physical modeling. We consider the case of a single crack with rough 
faces separating two media f2o and fix, which are both linearly elastic and isotropic, 
taking p to denote the density and c to denote the elastic speed of the compressional 
waves. These parameters are piecewise constant and may be discontinuous around 
the crack: (po, cq) if x S Oq, (pi, Ci) if x G Vt\. The media are subject to a constant 
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static stress p. At rest, the distance between the planes of average height is £o(f>) > 
(figure O left). 

Elastic compressional waves are emitted by a singular source of stress at x — x s < 
a in where a is a median plane of the actual flaw surface. The wave impacting 
a gives rise to reflected (in Oq) and transmitted (in Sli) compressional waves. These 
perturbations in SIq and fli are described by the 1-D elastodynamic equations 



dv 
~dt 



da 

dx' 



da 
~dt 



dv 



pc 



S(t)S(x - x s ), 



(2.1) 



where S(t) denotes the causal stress source, v — is the elastic velocity, u is the 
elastic displacement, and a is the elastic stress perturbation around p. The dynamic 
stresses induced by the elastic waves affect the thickness £ (t) of the crack (figure 12. 1[ 
right). The constraint 



£ = £o + M > f o - d > 



(2.2) 



must be satisfied, where [it] = u + — u~ is the difference between the elastic displace- 
ments on the two sides of the crack, and d(p) > is the maximum allowable closure 
@]. We also assume that the wavelengths are much larger than £, so that the propa- 
gation time across the crack is neglected, the latter being replaced by a zero-thickness 



interface at x 



[u] = [u(a, t)] = u(a + , t) — u{a , t). 
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Fig. 2.2. Sketch of the nonlinear relation between the stress and the jump of the elastic dis- 
placement at a. Left row: model 1 \2. 7] ), right row: model 2 12.86 . 

Two independent jump conditions at a need to be defined to obtain a well-posed 
problem. The discontinuity of a is proportional to the mass of the interstitial medium 
present between £Iq and Q± |24j . Since the crack is dry and contains only air, the 
density of which is much smaller than po or pi, the elastic stress is assumed to be 
continuous: 

[a(a, t)] = => a{a + , t) = a(a~ , t) = a*(t). (2.3) 

Establishing the second jump condition is a more complex task. Experimental and 
theoretical studies have shown that u is discontinuous, and that the discontinuity is 
proportional to the stress applied. The linear model has often been considered [21] : 

a*{t)=K [«(<*,*)], (2.4) 

where K is the interfacial stiffness. Welded conditions [u(a, t)] = are obtained if 
K — > +oo. However, the linear condition (|2.4p violates (12. 2\i under large compression 
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loadings conditions: cr*(i) < —K 5 =>■ £ < £o — 5. The linear condition (|2.4p is 
therefore realistic only with very small perturbations. With larger ones, a nonlinear 
jump condition is required. 

To develop this relation, it should be noted that compression loading increases 
the surface area of the contacting faces. A smaller stress is therefore needed to open 
than to close a crack; an infinite stress is even required to close the crack completely. 
In addition, the constraint ()2.2ll must be satisfied, and the model must comply with 
(|2.4[) in the case of small stresses. Lastly, concave stress-closure relations have been 
observed experimentally [17] . Dimensional analysis shows that the general relation 

a*{t) =KdF([u(a, t)]/d) (2.5) 

is suitable, where T is a smooth increasing concave function 

T : } - 1, +oo[^] - oo, Fxaaxi, um F(X) = -oo, < T max < +oo, 

X ^~] (2-6) 
= 0, T (0) = 1, T < < T . 

Two models illustrate the nonlinear relation (|2.5p . First, the so-called model 1 pre- 
sented in [2l [4] is 

Secondly, the so-called model 2 presented in [T7] is 

a*(t) =K dln(l + [u(a, t)]/d) T{X) = ln(l + X), F max = +oo. (2.8) 



These two models are sketched in figure 12.21 The straight line with a slope K tan- 
gential to the curves at the origin gives the linear jump conditions 



2.2. Numerical experiments. Here, we describe the influence of the nonlinear 
jump condition (|2.5p on the wave scattering. For this purpose, we consider a single 
crack described by model 1 (|2.7p . a sinusoidal source S with a frequency of 50 Hz, 
and parameters 

J Po = Pi = 1200 kg.nT 3 , K = 1.3 10 9 kg.m" 1 ^" 2 , 
\ c = ci = 2800 m.s- 1 , d = 6.1 10~ 6 m. 

The amplitude vq of the incident elastic velocity ranges from 10 -4 m/s to 5 10~ 3 m/s. 
This latter maximal amplitude corresponds to a maximal strain e = vq/cq ~ 10 -6 , 
so that the linear elastodynamic equations (|2.ip are always valid pQ. The linear 
first-order hyperbolic system (|2.ip and the jump conditions (|2.3p and (|2.7p are solved 
numerically on a (x, t) grid. For this purpose, a fourth-order finite-difference ADER 
scheme is combined with an immersed interface method to account for the jump 
conditions |15| . At each time step, numerical integration of v also gives u. 

Figure 12.31 shows snapshots of u after the transients have disappeared and the 
periodic regime has been reached. The mean values of the incident and reflected dis- 
placements (x < a) and the transmitted displacement (x > a) are given by horizontal 
dotted lines. With vq = 10~ 4 m/s (a), these mean values are continuous across a. 
At higher amplitudes, a positive jump from a~ to a + is observed. This jump, which 
amounts to a mean dilatation of the crack, also increases with vq (b,c,d). 
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Fig. 2.3. Model 1 \2. 7p : snapshots of the elastic displacement, for various amplitudes vo of the 
incident elastic velocity: 10~ 4 m/s (a), 10~ 3 m/s (b), 2 10 — 3 m/s (c) and 510 -3 m/s (d). The 
vertical solid line denotes the location a = 200 m of the crack. The red and green dotted horizontal 
lines denote the mean value of the elastic displacement on both sides of a. 



These properties can be more clearly seen in figure 12.41 where the numerically 
measured time history of [u] is shown, with «o = 2 10~ 3 m/s (a) and v o = 5 10 -3 m/s 
(b). It can also be seen from this figure that the maximum value of [u] increases with 
vq, whereas the minimum value of [u] decreases and is bounded by — d, as required by 
(|2.2p . These findings will be analysed whatever T in (|2.5p in the following sections. 



Distortion of the scattered fields can also be observed in figures 12.31 and 12.41 
increasing with the amplitude of the forcing levels. It is beyond the scope of this 
paper to sudy this classical nonlinear phenomenon. A local analysis in the case of 
model 1 (|2.7p was presented in [TOj . 



3. Model problem. 



3.1. Link with an ODE. To explain the numerical findings, we look for an 
evolution equation satisfied by [u(a, t)]. Considering the elastic wave emitted by the 
source S in (|2.ip . which impacts a single crack modeled by (|2.3|) and (|2.5p . leads to 
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Fig. 2.4. Time history of [u(a, t)], at two incident elastic velocities: do = 2 10 3 m/s (a) and 
vq = 5 10 — 3 m/s (b). The horizontal dotted line denotes —d. 



the following proposition. 

Proposition 3.1. The nondimensionalized jump of displacement y = [u(a, t)]/d 
satisfies the nonautonomous scalar differential equation: 



^7 = -F(y) + 1 ^^S(t/P), y(Q)= j/o = 0, 
dt pd p Cq 



(3.1) 



where (3 = i^((poCo) 1 + {pi ci) _1 ) > 0. The time t in hS.l}) has been successively 
shifted: t — (a — x s )/cq — > t, and rescaled: t — > fit. 

Proof. We adapt a procedure described in J22[ [5l [13] . The elastic displacement 
interacting with the crack is 



u(x, t) 



uj(t — x/cq) + un(t + x/cq) if x < a, 
Urit — x/cx) if x > a, 



(3.2) 



where uj, ur and ut are the incident, reflected and transmitted displacements, re- 
spectively. The elastic stress deduced from (|3.2p is Q] 



a(x, t) 



~p coUj(t - x/c ) + p a c u R (t + x/cq) if x < a, 



if x > a. 



(3.3) 



(3.4) 



—pi ci u T (t — x/ci) 

Calculations based on (|3.2p and (|3.3p yield 

a(a~, t) = -2 p c uj(t - a/c ) + p c u (aT, t), 
a(a + , t) — —pi ci u (a + , t). 

At a, the jump in the displacement \ an d the weighted displacement ip are introduced 

x(t) = [u(a, t)] = u(a + , t) — u(a~, t), 

iK*) = 7^ («(«". t) + ^u(a+, t) 
2 V Poc 
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or, alternatively, 



«(« > *) = ; W) + ; X(t), 

Po Co + Pi Cl p c + pi Ci 

/ - ,\ 2p c , . pici 
«0 , *) = V>(*) XW- 

PO Co + Pi Ci po C + pi Ci 



(3.5) 



The continuity of a (|2.3p means 



= %(t-a/c ). (3.6) 



The nonlinear relation (|2.5[) and the first equation of fl3.4p give 

/ / Kd 

u (a , t) = 2 %(t — a/co) H T{x/d). 

Po Co 

Using (|3.5[) and (|3.6p , we obtain after some operations 

= -/3dF(x/d) -2uj{t - a/co), 



where j3 is defined in proposition 13.11 Classical calculations of elastodynamics yield 
the incident elastic velocity generated by S and impacting the crack [T] 

v(x, t) = uj(t - x/cq) = - 1 2 S (t - - — — ) . (3.7) 

2poCq V c J 

The time shift t — (a — x s )/cq — > t therefore gives the differential equation 

d £ = -PdF{x/d) + -^S{t). 
dt p cfi 



Nondimensionalization y = x/d and time scaling t — ► f3t give the ODE (|3.1j) . Since 
the source S is causal and x s < a, then y(0) = 0. □ 

3.2. Sinusoidal forcing. From t = 0, the source is assumed to be sinusoidal 
with the angular frequency f2 

S(t) = 2v poc 2 sinfti, (3.8) 

which results in an incident elastic velocity with the amplitude vq (|3.7p . The source 
(|3.8p is injected in (|3.ip . Then, setting the nondimensionalized parameters 

. 2v H 2tt 27T/3 

A = 7TJ' w = 7?' T= — = ~FT' 

pa a u \l to a\ 

y = — 2 — ' J(y) = -F{y)> /»nn = -^w, 

the model problem becomes a nonautonomous differential equation with sinusoidal 
forcing: 

( dv 

= f(v) + A s ' mcjt = f)> ^ > o, o < a < +oo, 

/ :] - 1, +oo[-s-]/ min , +oo [, lim f(y) = +oo, -oo < / min < 0, 

(3.10) 

/(0) = 0, /'(0) = -l, /(y)<0</"(y), 
{ y{0)=y e] - 1, +oo[. 
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For the sake of generality, contrary to what occured in proposition 13.11 yg can differ 
here from 0. The properties of / in (|3.10p mean that the reciprocal function f^ 1 
satisfies 

rM/nun, +00H-l,+00[, km /- 1 (tf) = -l, 

v^+°° (3.11) 

r 1 (o) = o, (r 1 ) <o< (r 1 ) . 

In model 1 (12. 7|) . / is involutive and / m i n = — 1. The function / and its derivatives 
are 

f(y) = -Tzv,> f'w = ~ ?i 1 \2 ' f"( y ) = nl ^ /"fa) = - ttz ^4 - 

1+2/ (1 + 2/) (1 + 2/) (1 + 2/) 

In model 2 (|2.8[) . / m i n = — oo, and the function / and its derivatives are 

f(y) = -ln(l + y), f\y) = /'(„) = —i-j, /"(„) = -^-^3. 

1 + 2/ (1 + 2/) (1 + 2/) 

In the next sections, solutions of (|3.10[) and related ODE will be presented. These 
solutions are computed numerically using a fourth-order Runge-Kutta method, with 
model 1 (|2.7p . The parameters are those used in section [2~2l with nondimensionaliza- 
tion (|3 . 9(1 . In addition and for the sake of clarity, the dependence of solutions on the 
parameters A and uj is omitted, except when necessary. Lastly, we use an overline to 
denote the mean value of a function during one period [0, T]\ given s(t), 

1 f T 

4. Preliminary results. 

4.1. Existence and uniqueness of a periodic solution Y. The T-periodic 
isocline of zero slope deduced from (|3.10| is 

I%(t) = TH- A sin LJt). (4.1) 

PROPOSITION 4.1. There is a unique T-periodic solution Y(t) of S3.10\) . This 
solution is asymptotically stable. 

Proof. For the sequel, we denote Y = Y(0) the unique initial data such that the 
solution of (|3.10p is T-periodic. Three cases can be distinguished. 

Case 1 : A = 0. Equation (13.10|) becomes a scalar autonomous equation with the 
steady state solution Y = 0. Since / is a C 1 function and / (0) = — 1, the fixed point 
is asymptotically stable. In the case of model 1 (|2.7[) . the exact solution is a known 
special function: ye v — yoe Va ~ t =>■ y(t) — C (j/o e 1 " 1- *), where C is the Lambert 
function [6]. 

Case 2 : < A < |/ m in|- The T-periodic isocline (|4.1|) is defined on M., with 
lower and upper bounds f~ 1 {A) < and / _1 (— A) > (figure l4TTT -a). Below Iq, we 
obtain 4| > 0; and above Iq , we obtain -j| < 0. Consequently, the null slope of the 
horizontal line f _1 (A) is smaller than F(t, y) when these curves intersect, except at 
t = 3T/4, where they are both equal to zero: f^ 1 (A) is called a weak lower fence 
for ()3.10|) [TO] . Similarly, F(t, y) < (/ _1 (— A)) = 0, where equality occurs only at 
T/4: / _1 (— A) is called a weaA; upper fence for (|3.10|) . Since j/ 1— > i^(t, y) is Lipschitz 
continuous, /~ 1 (±A) are nonporous fences: once a solution has crossed these fences, 
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it cannot sneak through them [10] . Based on the funnel's theorem, the subset of the 
(t, y) plane defined by 

K = R+xK Q , K =[f- 1 (A) > r x (-A)], 

is a funnel [10] : once a solution has entered K, it stays inside. Defining the subsets 

U+ = {t> 0; -1< y < r\A)} , U~ = {t > 0; y > f~ 1 (-A)} , 

the solution of (|3 . 1 0[) increases in U + and decreases in U~ as t increases. In addition, 
zero slope 4f = in occurs only at T/4 and 3T/4. The solution therefore 



enters K in finite time, and Kq is invariant under the flow tp : 



1, +oo[ 



] — 1, +oo [; (t, yo) I— > y(t) = ip(t, y ). The solutions of (|3.10p are therefore always 
bounded, which ensures that the solutions are global solutions in time. 

The invariance of Kq ensures IL(Ko) C Kq, with the Poincare map LT :] — 1, +oo[^ 
] — 1, +oo[; yo i > n(yo) = f(T, yo)- Since II is continuous on K, it has at least one 
fixed point Y . The solutions Y(t) = (p(t, Y ) and Y T (t) = tp(t + T, Y ) both satisfy 
the same ODE with the same initial value: Y T (0) = <p(T, Y ) = IL(Y ) = Y a = Y(0). 
Therefore Y(t + T) = Y(t) at all values of t, which proves that Y(t) is T-periodic. 

(a) , 

upper fence f " [ (-A) 





Fig. 4.1. Phase portrait of (3.10}) with model 1 \2. 7] ), at A = 0.6 (a) and A = 1.2 (b): isocline 
of zero slope Iq and periodic solution Y. Dashed horizontal lines: lower fence / _1 (j4) < 0, upper 
fence A) > (a). The dotted horizontal line stands for y = —1. In (b), the dotted vertical 

lines stand for asymptotes at T\ and T2 . 

Let us consider another T-periodic solution Y(t) of (|3.10|) . with Y(0) = Yq. As- 
suming that Y > Y yields Y{t) > Y{t) at all values of t <E [0, T], hence f(Y(t)) < 
f(Y(t)). The integral forms of Y(T) and Y(T), together with the T-periodicity, give 

Y T -Y T = Y0-Y0+ (f(Y(t)) - f(Y(tj) dt<Y - Y , 
= Y -Y . 

This results in Yq — Yo < Yq — Yo, which is impossible. The opposite case Yo < Yq can 
be handled in the same way, which proves the uniqueness of the T-periodic solution. 
Lastly, the properties of / in (|3.10p mean that 

a °= / jr(*><p(t>Vo))dt= f f\<p(t,ya))dt<0, (4.3) 
Jo d y Jo 
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therefore Y(t) is asymptotically stable (see e.g. theorem 4-22 in [5]). 

Case 3 : A > |/ mm |. This case can occur only if / min > — oo. The isocline (|4.ip is 
defined outside the two vertical asymptotes at T\ and Ti 

0<T 1 = iarc8in(fei) <|<T 2 = |-T 1 . (4.4) 

(figure l4~Tl b). The set [f^ 1 (A), +oo[ is invariant under the flow but not compact. 
To prove the existence of solutions at all time, let us examine the subset of the (t, y) 
plane where the solution can blow up: T = {0 < t < T/2; y > 0}. In T, y > and 
f(y) < 0, hence 4| < A smut. By integrating, we obtain 

f T/2 A , 2A 
< y < yo + A sin cut at = yo H < +oo, 

Jo ^ 

which proves that the solution is global in time. 

To obtain an invariant set, we take y > 0, which means that /min < f(y) < and 
thus / m i n + A smut < 4| < A sinuit. As long as y > 0, integration will ensure that 
y is bounded by a lower solution and an upper solution 

A A 

yo + /min t H (1 - coswi) < y(t) < y + — (1 - cosc^) . 

uj u> 

Taking yo > |/min| T, for instance yo = 2 |/min| T, ensures that y > for all values of 
t E [0, T) and that y(T) < yo- Therefore, upon defining the subset of the (t, y) plane 

K = M+ X K , K = [/ _1 (A), 2 |/min| T] , 

then i^o is invariant under the Poincare map, yielding II (/<o) C Kq. The same argu- 
ments as those given in case 2 hold here, which proves the existence, uniqueness and 
asymptotic stability of the T-periodic solution. □ 

As seen in the previous proof, the phase portrait of (|3 . 1 0|) depends on A. However, 
the evolution of Y with t is the same at all values of A, as seen in figure 14.11 and 
proved in section [5] Some auxiliary solutions are first introduced. These solutions 
make it possible to investigate how the attractive periodic solution Y evolves with the 
parameter A. 



4.2. Auxiliary solutions. The first derivative of Y with respect to A is intro- 

dY 
OA' 



duced: Z(t, A) — ^j. Applying the chain-rule to ()3.10j) shows that the T-periodic 



solution Z satisfies 



*<M)-*-g. 

The isocline of zero slope of (|4.5|) is deduced from (|3.10|) : 



^ = f\Y)Z + smut, 

dt ... (4.5) 



r y , . smut smut ., „. 

7 o (*) = - ~7~~777~7 = — ■ (4-6) 



/ (Y(t)) 



f (Y(t)) 



Below 1$ , we obtain > 0; and above Iq , we obtain < 0. The phase portrait 
of (|4.5p obtained in the case of model 1 is displayed in figure 14.21 which shows the 
qualitative properties of Z that will be described in the next lemma and proposition. 
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(a) , 

isocline of zero slope 




isocline of zero slope 




Fig. 4.2. Phase portrait of {4.5}) with model 1 {2.7}) , at A = 0.6 (a) and A = 1.2 (b): isocline 
of zero slope Iq and periodic solution Z. 



Lemma 4.2. Setting 



the initial value of the T '-periodic solution Z in \4-5\ ) satisfie 



B(s) = - / / (Y(r))dr>0, 



Zn = 



e B ^ 1 Jo 



e B ^ sin uisds < 0. 



(4.7) 



Proof. Given Y(t), the ODE (|4. 5|) is linear, which leads to the closed-form ex- 
pression 



Z{t) = + e B{s) sinus ds^j e 



B(t) 



T-periodicity of Z gives D and proves the equality established in (|4.7p . Since s i— > e s ( s ) 
is a strictly positive and increasing function of s, we obtain the following bound 

fT fT/2 pT 

/ e B{s) sinus ds < e B{T/2) I sinus ds + e B{T/2) \ sinusds = 0. 

JO JO Jt/2 



Substituting this inequality into (|4.7|) proves that Zq < 



□ 



Proposition 4.3. In [0, T], the T -periodic solution Z of |^.5[ j /ias two roofs tz x 
and tz 2 , which are ordered as follows: 

0<t Zl < T/2 < t Z2 < T. 

Therefore Z < in [0, tz 1 [L)]tz 2 T], and Z > on ]tz 1 , tz 2 {- In the limit of null- 
forcing case A = 0, Z and its roots are determined analytically: 



Z{t, 0) = 



VT+ur 



sva(u)t — 6), 9 = arctanw, 



Wo) 



6 T 
u> + 2' 



(4.8) 
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Proof. T-periodicity of Z and integration of (|4. 5|) over one period gives 

/ /(Y(r))Z(r) = 0. 
Jo 

Since / < 0, the sign of Z changes in [0, T]. From lemma [4T2l and the phase portrait 
of (|4.5j) . it can be deduced that Z must cross the t-axis twice, at tz t and tz 2 in 
]0, T[. Let us assume tz 1 > T/2, where the isocline (|4.6[) is negative (figure ■ 
Then Z will intersect Jq twice at negative values and will never cross the i-axis, 
which is impossible. After tz 1 , Z increases and intersects Iq with a zero slope at 
tzi < t < T/2. It then decreases but cannot intersect again Iq as long as I§ > i.e. 
at t < T/2. This means that t Z2 > T/2. 

Taking A = results in Y = 0, as stated in the proof of proposition 14. II fcase 1), 
and hence (|4.5[) becomes the linear ODE 

Z = —Z + sinujt, 
the T-periodic solution of which is given in (|4.8|l . □ 




Fig. 4.3. Phase portrait of uraift modeZ J p.7| ), at ^4 = 0.6 (a) and A = 1.2 ('6J: isocline 

of zero slope 1^ and periodic solution W. 



The second derivative of Y with respect to A is also introduced: W(t, A) — '' 1 



Applying the chain-rule to (|4.5[) shows that the T-periodic solution W satishes 



dA 1 



^-=f"(Y)Z 2 + f'(Y)W, 
W(0, A) = W = 



(4.9) 



The isocline of zero slope of W, which is deduced from (|3.10p and (|4.9j) . is 



/to) 



Z 2 (£) > 0. 



(4.10) 
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This isocline, which is shown in figure 14.31 is equal to zero only at tz t and tz 2 , when 
Z vanishes ( proposition ^. 3|) . Below I™ , we obtain ^j- > 0; and above , we obtain 

dt V u ' 

Proposition 4.4. The T -periodic solution W of j^.ffp is strictly positive: for all 
t e [0, T], W(t, A) > 0. 

Proof. If W < 0, the phase portrait of (|4~9| entails that W{T) > Wo, which is 
contradictory with the T-periodicity of W: hence, W > 0. We then consider 

r w=f'(Y) w , 

\ < w(0) < W . 

Since /" > and Wo > w(0), w is a lower solution of W, hence W(t) > w(t) at all t. 
However, the T-periodic lower solution is w = 0, which completes the proof. □ 

Lemma 4.5. In the limit case of null forcing, the T-periodic solution W of h4-9ty 
contains only the harmonics sin 2tot and cos 2ujt, and has a non-null mean value 

no) = ^j^. (4.1D 

Proof. We inject Y = into (|3T0|) and (|49|) . Using ([48)1 gives 

W'(t, 0) = -^-^ sin 2 (ujt-e) -W(t, 0). 
1 + OJ 

Integrating the latter equation over [0, T] gives the result. □ 



Lastly, the third derivative of Y with respect to A is introduced: X(t, A) = f-^r- 
Applying the chain-rule to (|4.9[) shows that the T-periodic solution X satisfies 

^ = /'" (Y) Z 3 + 3 /" (Y) Z W + /' (Y) X, 
'" 0* Yo ( 4 - 12 ) 



X(0, A)=Xo = 



dA 3 



Lemma 4.6. In the limit case of null forcing, the mean value of the T-periodic 
solution X of is X(0) = 0. 

Proof. Upon substituting Y = into (|3.10|) . then (14.121) becomes 

X'(t, 0) = f"'(0)Z 3 (t, 0) + 3f"(0)Z(t, 0)W(t, 0)-X(t, 0). 
Integrating this differential equation over [0, T] and the T-periodicity of X give 



X(0) = f (0) Z(t, 0) 3 + 3 / (0) Z(t, 0) W(t, 0). 



In (|4.8|) , Z(t, 0) is a sinusoidal function with a null mean value, hence Z(t, 0) 3 = 0. 
Likewise, (|4.8p and the properties of W(t, 0) stated in lemmag3]give Z(t, 0) IV (t, 0) = 
0, which proves the result obtained. □ 
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0.5 1 1.5 2 2.5 3 

A 



Fig. 5.1. Evolution of Yq with the forcing amplitude A in the case of model 1 \2. 7| ). See 
proposition 15. 1\ 



5. Qualitative properties of Y. 

5.1. Phase portrait of Y. The properties of Z and W provide qualitative 
insight about Y, stated in the next propositions and corollary. 

Proposition 5.1. The initial value Yq in Ii3.10\) is a strictly decreasing convex 
function of A, which satisfies 



-1 < max ( f~ l (A), 



A. 



< Y {A) < 0. 



Equality Yq = occurs only at A — 0. At low forcing levels, we obtain the local 
estimate 



Y 



A. 



1+UJ Z 



0{A 2 ). 



Proof. Propositions and SU yield £f = z{0, A) < and ^ = W(0, A) > 
0, which proves the convex decreasing of Yq. Since 1q(0) = Y(0, 0) = 0, then Yq(A) < 
0, which gives the upper bound of Yq. 

We now prove the lower bounds of Yq. In the proof of proposition 14. li it was 
established that f~ 1 (A) is a lower nonporous fence, which means that Y(t) > f~ 1 (A), 
and hence that Yq > f^ 1 (A). On the other hand, lemma l4~2l and proposition 14.41 
show that A i— » Yq(A) is a convex function, which reaches a maximum at A = 0, with 
Y"o(0) = 0. Yq(A) is therefore always above the straight line going through the origin 
with slope A f^-(O): Y > AZ (0), where Z (Q) is deduced from 

At A <C 1, a Taylor expansion of F is written as follows 

dYn 



Y a (A)=Yo(0) + A^(0) + 0(A 2 ). 
E A = 0, then F = Y = 0, and gl]) again provides Z o (0) = |^(0) 







The evolution of Yq with A is shown in figure 15.11 Critical times in Y are defined 
and estimated in the following proposition and corollary. 

PROPOSITION 5.2. In [0, T], the T -periodic solution Y of \3.1U\) , which has two 
roots ty t , ty 2 , is maximum and positive at ty max , and minimum and negative at ty min ■ 
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Fig. 5.2. Critical times in Y and Z, in the case of model 1 \2. 7] ). Snapshot at A = 0.6 (a): the 
vertical dotted lines denote successively tz 1 , t Yl , iy max > W 2 > tz 2 an d ty ■ ■ Parametric study with 
A (b): a horizontal dotted line has been added at T/2, and four horizontal dashed lines have been 
added at 9/uj + iT/i (i = 0, ...,3) to show the locations of critical times for A — > 0. See proposition 
\5.2\ and corollarv \5.3\ 

The critical points are ordered as follows: 

< t Yl <T/2< ty 2 < T, t Yl < iy max < iy 2 < t Yoxin , 
T/4 < <y max < T/2, 3T/4 < t Yxaia < T. 

Wit 4cl and using 9 defined in J^.£[ ), we also obtain the following local estimates 

Y(t, A) = J— sin (ut -6)+0 (A 2 ) , 

, i . 9 . i . 9 T , i. T . , , 9 3T 

t Yl (Q+) = - t Y _(0+) = - + -, t Y2 (0+) = - + -, < Ymin 0+ =- + — . 

u> u> 4 u> 2 uj 4 

(5.1) 

Proof. The phase portrait of (|3 . 1 0|) and the proposition 15.11 mean that Y(t) must 
cross the t-axis twice, at t Yl and ty 2 in ]0, T[. At t > T/2, the isocline Iq (|4.ip 
is negative (figure |4~T|) . Consequently, t Yl > T/2 means that Y crosses Iq twice at 
negative values of Y, preventing Y from crossing the t-axis, which is impossible. 

The slope of Y is null at t Ytaax , which means that Y crosses Iq. This occurs 
only when the slope of Iq is negative, between T/4 and 3T/4. Since 1^ must also 
be positive at the intersection point, it is possible only when T/4 < iy max < T/2. Y 
subsequently decreases but cannot cross the t-axis as long as Iq positive, which means 
that iy 2 > T/2. 

Lastly, the slope of Y is equal to zero at ty min , which means that Y crosses Iq . 
This is possible only when the slope of Iq is positive. Since Iq must also be negative, 
this means that £y min > 3 T/4. 

To prove the local estimates, a Taylor expansion of Y is written as follows 

Y(t,A) = Y(t,0)+A^(t,0) + O(A 2 ), 
= AZ(t,0)+O(A 2 ). 
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Equation (|4. 8|) and straightforward calculations complete the proof. □ 

COROLLARY 5.3. The following properties hold for all A > 0: 
(i) [t Yl , ty 2 ] C \tzx-, tz 2 ]> or equivalently Z(t Yl , A) > and Z(ty 2 , A) > 0; 

(11) ^X <0 ' -dA >0; 

e t e t 

H t Yl < - < -, ty 2 >- + 77; 

T 

(iv) if / m in < oo, i/ien lim = — . 

A^ + oo Z 

Proof. 

Property (i). Two subsets of the i-plane are defined: 

1(A) = {t e]0, T[; F(t, A) > 0} , J(A) = {f e]0, T[; Z(t, A) > 0} . 

Equations P~g)) and (|5TTj) prove that 1(0+) = J(0). Phase portraits of ([3~TU)l and 
(|4.5p show that / and J are open intervals: 

I{A)=]t Yl ,ty 2 [, J(A)=]t Zl ,t Z2 [. 

Let us consider two forcing parameters A\ and A 2l with A2 > Ai > 0. Proposition 
gives |f > 0, and hence Ai) < A 2 ). Taking t G J(Ai) gives < 
Z(t, Ai) < Z(t, A 2 ), and therefore t S J(^2): the interval J increases strictly with 
A. To prove 1(A) c J(^), we take £ e ^(-A), and hence Y(t, A) > 0. Two cases can 
be distinguished: 

• t (£ 1(0+) Y(t, 0+) < 0. The formula 

f A 

Y(t, A) = Y(t, 0+) + / Z(t, s) ds 



Jo+ 

and the monotonic increase in J with A yields Z(t, A) > 0, hence t £ J(A); 
• t € ^(0 + ). Arguments mentioned above yield £ G J(0) C J (A), which com- 
pletes the proof of (i). 
Property (ii). First we differentiate Y(t Yl , A) = in terms of A: 

dY dt Yl 9Y 
lftll~A + ~dA~ ■ 

The time derivative is replaced via (|3.10p 

(f(Y(t Yl , A)) + A sin ujt Yl ) ^ + Z(t Yl , A) = 0. 
Since f(0) = 0, we obtain 

dt Yl _ Z(t Yl , A) 
dA A sin ujt Yl 

Property (i) and proposition 15 . 21 then prove -pj- < 0. A similar derivation gives 

dt Y2 _ Z(t Y27 A) 
dA A sinujt Y2 
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Property (ii) and proposition 15.21 give Z(ty 2 , A) > and smujty 2 < 0, and hence 

OA ^ U ' 

Property (iii). The bounds of ty 1 and ty 2 follow from (|5 . 1 [) and property (ii). 
From the definition of 9 (I4.8[) and the property of arctan, we also obtain 

9 1 7T _ T 

u! ~ u 2 ~~ T ' 

which provides an additional bound for ty 1 ■ 

Property (iv). If A > |/ m in|, Y crosses the isocline Iq at iy max > T 2 : see (I4.ip . 

(|4.4j) and figure 14. U(b). The bounds of proposition I5T21 ensure T 2 < iy max < T/2. 
Using lim T 2 = T/2 completes the proof. □ 

A— > + oc 

Numerical computations illustrating proposition 15.21 and corollary 15.31 are shown 
in figure [ST2J No theoretical results have been obtained about - and d g^g 1 " . The 
numerical simulations indicate that ty max strictly increases with A, whereas iy min first 
decreases up to a critical value, and then increases (b). Numerical simulations also 
indicate that tz 2 tends asymptotically towards £y min at large values of A. 




Fig. 5.3. Evolution ofY with the forcing amplitude A, at < A < 2 (a) and < A < 20 (b), 
in the case of model 1 \2. 7] ). See theorem \5.4\ 

5.2. Mean dilatation of the crack. As stated in the introduction, the main 
aim of this study was to prove that a jump in the mean elastic displacement occurs 
around the crack, as observed numerically. This jump amounts to a mean dilatation 
of the crack. The next theorem addresses this typically nonlinear phenomenon. 

Theorem 5.4. The mean value of the T '-periodic solution Y in \3.1U\) is positive 
and increases strictly with the forcing amplitude: 

- dY 

r>o, ->o. 

At small forcing levels, the following local estimate holds: 

Y=LM T ^_ +0 (A i ). (5.2) 
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Proof. T-periodicity of Y and sinwi in (|3.10p yields f(Y) = 0. Applying Jensen's 
inequality to convex / gives f(Y) < f(Y) = 0. Since /(0) = and / < 0, then 
Y > 0. PropositionHUand g^D yield Z(0) = and W = §| > 0, and hence Z(A) = 



_ OA 

§^(A) > 0, which proves the second inequality. At A <C 1, a Taylor expansion gives 



«v /2 a2y 4 3 r) 3 Y 

= Y(0) + AZ(0) + ^-W{0) + ^-X(0) + O(A 4 ). 
2 b 

From ([5TT]) and (|4~8|) . we obtain F(0) = and Z(0) = 0. Lemmas 1431 and l4~6l then 
give the local estimate of Y. □ 



The evolution of Y with A is presented in figure 15.31 With A = 2, the relative 
error between Y and its local estimate (|5.2[) is less than 5%. Figure IQl -fa) may be 
rather misleading as far as moderate values of A are concerned: it might seem to 
suggest that Y is always greater than the local estimate in (|5.2p . This is not in fact 
the case with greater values of A: the position of Y relative to its local estimate is 
not constant (b). 

5.3. Maximum aperture of the crack. Let 

^max = Y(tY w ) = SUp Y {t) 

te[o,T] 

be the maximum value of the T-periodic solution Y . Wc introduce a technical lemma. 

Lemma 5.5. Let h £ C 3 (I x A, R), where I and A are open subsets o/R. The 
following properties are assumed to hold whatever the value of A in A: 

(i) h + (X) = sup/i(z, A) is reached at a single point z + (A) G I; 

zei 

(ii) ^4(z + (\), A) < 0; 
o z 

T 9h n d 2 h n 

(m) yzeI '8X >0 '0>? >0 - 
Then h+(X) e C 2 (I, R), and A h-> /i+(A) is a sinc% increasing and convex function. 

Proof. The definition of z + in (i) gives 



|^(z+(A), A) =0. (5.3) 



In addition, (ii) ensures that ^§(z + (A), A) ^= 0. The implicit function theorem can 
be applied to §f e C 2 {I x A, R): therefore, A ^ z+{\) e C* 2 (A, R), and h+(A) = 

h(z+, A) e c 2 (/"r). 

Based on the first property in (iii) and (|5.3[) . 

dh + dh, , dh, , .-,dz + 

which proves that fo + is a strictly increasing function of A. On the other hand, 
differentiating (|5.3p in terms of A gives 

■(*+, A) + — A) — = 0. (5.4) 



d\d z dz 2 dX 
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Using (ii), the second property in (hi), and (|5.4p . we obtain 
d 2 h+ 



d\< 



-(A) 



d 2 h . , d 2 h , , ^.dz + 



9A 
d 2 h 



dzd A 



dA ' 



which proves that h + is a strictly convex function of A. □ 



Theorem 5.6. Y max is a strictly increasing convex function of A, and the fol- 
lowing bounds hold for all A: 



A A 1 T A 2A 

r 7; i r. o o /mill I < max < 



(5.5) 



The first lower bound of Y max is also a local estimate for small forcing levels A <C 1 

A 



vT 



0(A 2 



(5.6) 



5 

t 4 

> 

E 

§ 2 



(a) 


upper bound 






Ymax/ 




local estimate 



0.4 0.8 1.2 1.6 

A 




Fig. 5.4. Evolution of y ma x with the forcing amplitude A, at < A < 2 ('aj and < A < 20 
(b), in the case of model 1 \2. 7| ). See theorem 1 5. 61 



Proof. To analyse the evolution of F max in terms of A, we proceed in three steps. 
First, the phase portrait of Y and the proposition 15.21 show that Y(t, A) reaches its 
maximum value at a single time between T/4 and T/2. Secondly, the definition of 
ty mex an d proposition give 



d 2 y 



f (Y m 



dY 



dt 

Au coswiy- max < 



(£y max , A) + Aui COS LOty m 



Thirdly, corollary EH and proposition S3 yield > and > on [T/4, T/2]. 



The three hypotheses in lemma 15.51 are therefore satisfied, which proves the convex 
increasing of Y, nax with A. 
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The bounds of y max are proved by building upper and lower solutions of Y , as 
in the proof of proposition 14. 11 case 3. In ]ty 1 , iy max [, proposition 15.21 ensures Y > 0, 
which means that /,„;„ < f(Y) < 0, and hence / m ; n + A sinuit < 4y < A s'mcut. 
Integrating this inequality from ty 1 to iy maJI gives 

A A 

— (cOSWty — COS U)ty „) + {ty 1 ~ ty^J |/min| < Y ma x < — (COS Ulty 1 ~ COS U)ty ) . 
LU LU 

Proposition 15.21 and corollary 15.31 - (iii) give the bounds T/4 < iy max < T/2 and < 
ty 1 < 9/u, and hence 

{COS 9 < COS LOty 1 < 1 , 
< — COS U>ty max < 1. 

These inequalities, together with the definition of 9 in (|4.8[) . yield 
A 1 A , s 2A 

^=^= < — (cOSWty. — COS Ujty ) < . 

W VI + CJ 2 W " W 

Proposition l5 . 21 also gives ty — iy max > —T/2, which gives the upper bound and a first 
lower bound in (|5.5p . The latter lower bound is not always an optimum bound: when 
A < -j | /mini w \/l + w 2 , it is negative, whereas proposition 15.21 states that Ymax > 
at all values of A. We therefore take advantage of the convex increasing of Y max with 
A, which was previously proved. Since Y max (0) = 0, then Ym ax (^4) is always above the 
straight line going through the origin with the slope A Z(ty la ^(0 + ), 0). Proposition 
15.21 and equation (|4.8[) therefore give a second lower bound in (|5.5p . A second-order 
Taylor expansion of Y max (j4) proves that this lower bound is also the local estimate 
of Ymax when A < 1. □ 

The evolution of Ynax stated in theorem 15.61 is shown in figure [5^1 If |/ m i n | < oo 
and the two lower bounds given in (|5.5p intersect at A = -j |/ m in| jz^j Vl + 

At higher values of A, the second lower bound is more accurate than the first one. 




Fig. 5.5. Evolution of y m i n with the forcing amplitude A, with < A < 2 (a,) and < A < 20 
('oj, in t/ie case of model 1 \2. 7}). See theorem\5. 



5.4. Maximum closure of the crack. Let 

Ymin = Y(ty mm ) = mf Y(t) 
t£[0,T] 
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be the minimum value of the T-periodic solution Y. 

Theorem 5.7. The following properties are satisfied by Ymin-' 

(i) for all A, we obtain the bounds: 

-l<f-\A) <y min <0; (5.7) 

(ii) at small forcing levels, we also obtain the local estimate 

Y min = -^JL= + 0(A 2 ); (5.8) 
vi + ^ 

(iii) at infinite forcing levels, the crack is completely closed: 

lim Y min = -1. 

Proof. 

Property (i) . Proposition 15.21 shows that Ymi n is negative. As in the proof of 
proposition [STTJ the existence of the lower nonporous fence / _1 (A) means that Y mnl > 

f-HA). 

Property (ii). As in the proofs of proposition 15.11 and theorem 15.61 the local 
estimate when A -C 1 follows from a second-order Taylor expansion of Y and from 



Property (iii). First we consider the property f(Y) — used in the proof of 



theorem 15.41 Evolution of / in (|3. 10[) and the definition of t Yl , ty 2 in proposition 15.21 
give 



T 



hence 



= / f(Y(r))dT+ / f(Y(r))d T + / f{Y{r))dr, 

Jty x Jt Y2 

Yl \f(Y(r))\dr+ f \f(Y(r))\dr) + \f(Y(r))\dr, 
Jt Y2 J Jt Yl 



^ \f(Y(r))\dr+ ( T \f(Y(r))\dr) = F' \f(Y(r))\dr. (5.9) 



Y 2 / Y l 



I J 

Secondly, let us consider the ODE 

ds 

—— = —s + A smujt, 

dt ' (5.10) 

s(t Yl ) = 0. 

Straightforward calculations and property (iii) of corollarv l5.3l give 

s{t) = - A (\ sm(tot Yl - 0)\ e- {t - tY ^ + sin(^ - 6)) , (5.11) 

with 9 (|4.8p . Comparison between (|3.10p and (|5. 10[) shows that s is a lower solution 
of Y when Y > 0: Vi e [t Yl , ty 2 ], Y(t) > s(t). Since £ Fmax < ty- 2 , it follows from (EU) 
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and (j5TTTj) that 



J > / |/(Y(r))|dr> / |/(«(r))|dr, 



Y 1 •> t Yl 



> 



A 



VT+Z? Jt 

\/l+UJ 2 



(| sin(cjt Yl - 6»)| er {T - tY ^ + sin(cjr - 0)) dr, 



with 
A = 



| sm(ujt Yl -6)\ (l- e ( * y — "* y i^ + - (cos^ty-j - 0) - cos(wty max - 0)) . 
Proposition 15 . 21 and property (iii) of corollarv l5.3l give 

1 - e T ' 2 < 1 - e<**»~ < 1 - e 274-fl/« 
< | sin(wtyi - 0)| < sin0, 

2 cos0 1 . 1 — sin0 

< — (C0S(tL'ty 1 — 0) — COS(w£y max — 0)) < 



hence A is bounded independently of A. Finally, we obtain 

2A 



(! + «") 



J > 

J blows up when A — > +00. Equation (|5.9p means that I behaves in a similar way. 
Since Y < in ]0, ty x [U]ty a , T[ and given / in p.lOp . Y - must tend towards -1 when 
A — > +00, to make / blow up, which concludes the proof. □ 

The evolution of Ymm stated in theorem 15.71 is shown in figure 15.51 At moderate 
forcing levels (a) , one cannot distinguish between Ymi n and its lower bound. At higher 
forcing levels (b), one observes that Ymin is above its lower bound, as stated in theorem 
15.71 Although no rigorous proof has been obtained so far, the numerical simulations 
indicate that Y m i n decreases strictly as A increases. 

6. Generalization. 

6.1. Periodic forcing. Most of the results obtained in sections S] and [5] were 
based on the simple analytical expression of sinusoidal forcing. In the case of more 
general forcing, the key properties stated in proposition 15. 21 and corollary 15 . 31 are lost, 
which makes it impossible to obtain estimates such as theorems l5.6l and l5.7l However, 
some properties are maintained if two assumptions are made about S in ()3. 1|) : Xb- 
periodicity, and a null mean value. The latter assumption is physically meaningfull: 
as deduced from (|3.7[) . a non-null mean value of S results in an amplitude of uj that 
increases linearly with t. 

Since the techniques required are the same as those used in the above sections, the 
derivations will be shortened. Setting the parameters and the Fourier decomposition 
of the source as follows 

*'o = ^ max \S(t/0)\, A=^, y = M^M, f(y) = -F(y), 
2/9 Co te[o,T] pd d 

1 °° 2 7T 

9(t) = - o — \S(t/(3)\=} (a n sinnu;t + b n cosnujt), T = j3T , uj=—, 

2 Po c o v o ^ T 

(6.1) 
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we obtain the model problem 

' i% = m + Ag(t)=F(t,y), 

f :] - 1, +oo[— ^/min, +00 [, lim f(y) = +oo, -oo < / min < 0, 

< /(0)=0, /'(0) = -l, /(y)<0</"(y), ( 6 -2) 

fl(t + T) =.g(<), 5 = 0, |s(t)|<l, max Ut)| = 1, < A < +oo, 

te[o,r] 

k 2/(0) = yo S] - 1, +oo[. 
The isocline of zero slope of (|6.2p is 

7 a (t)=/- 1 (-A 5 (t)). (6.3) 

6.2. Periodic solution. The following result generalizes the proposition 14.11 
since the same notations and processes are used, the proof is only sketched here. 

Proposition 6.1. There is a unique T -periodic solution Y{t) of \6.2\) . This 
solution is asymptotically stable. 

Proof. Three cases can be distinguished. 

Case 1: A — 0. This case is identical to case 1 in proposition ^. II 

Case 2: < A < |/ m i n |. Based on the funnel's theorem, the compact set 

K =[r 1 (A), f~\-A)] 

is invariant under the flow of (|6.2p . and hence H(Ko) C Kq. A fixed point argument 
shows the existence of a T-periodic solution Y. Equations (|4.2p and (|4.3[) still hold: 
the uniqueness and attractivity of Y are straightforward consequences of / (y) < 0. 

Case 3: A > |/ m i n |. The isocline (|6.3p is not defined when Ag(t) > — / m i n . 
In these cases, the invariant set deduced from the isocline is the non-compact set 
[f^ 1 (A) 1 +oo [, which makes it impossible to apply the fixed point argument. To build 
a compact invariant set, we take y > 0, which means that fmin + Ag(t) < < Ag(t). 
As long as y > 0, integration gives lower and upper solutions of y 

VO + /mint +/ Ag(0d£ <y(t) <y Q + f Ag(t)d?. 
Jo Jo 

From the T-periodicity of g, it follows that 

Uo + fminT < y(T) < y Q . 

It is then only necessary to bound yo to establish that y > on [0, T]. Setting 

#min = min g{t) < 0, 
fe[o,T] 

the lower solution of y gives 

yo - 1 (|/ mi „| + A\ 9min \) < y (t), -it e [0, T]. 

Taking y > T(|/ m i n | + ^4|g m i n |) therefore gives y > 0. The compact set 

K = [rV), 2T(|/ min |+v4| 5min |)] 
is therefore invariant under the Poincare map, which completes the proof. □ 
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6.3. Auxiliary solution. Like in section^] the first derivative of Y with respect 
to A is introduced: Z(t, A) = ^-j. Applying the chain-rule to (|6.2p shows that the 
T-periodic solution Z satisfies 

^L = f'(Y)Z + g(t). (6.4) 
Lemma 6.2. In the limit of null forcing A — ? the mean values of Z and Z 2 are 

Z(0) = 0, ^2( ) = -^-^±n , (6 .5) 

2 n= i l + {nu) 



where a n and b n are the Fourier coefficients of the source k6.1\) . 
Proof. Null forcing gives Y = 0, hence (16 .4[) becomes 

d Z „ 

This equation is integrated over [0, T]. T-periodicity of Z and g = yield Z(0) 
The solution Z is therefore sought as a Fourier series with null mean value 

oo 

Z(t, 0) = (A n sinntut + B n cosnuit) . 

n=l 

Injecting this series in the ODE satisfied by Z(t, 0) and using (|6.ip provides 
a n + nub n b n -nuja n 2 2 + b 2 n 

A n = ; -T5-, B n = ; -tt , With A„ + B.„ = 



Now, Parseval's formula leads to the last equality in (|6.5p . □ 

6.4. Mean dilatation of the crack. The next result extends the theorem 15.41 
It shows that a positive jump in the mean elastic displacement still occurs. 

Theorem 6.3. The mean value of the T-periodic solution Y in \6. 2\) is positive 
and increases strictly with the forcing amplitude: 

— dY 
Y>0, ->0. 

At small forcing levels, the following local estimate holds 

F = /"(°)(i) ETT7 t \ + °(^ 3 )' ( 6 - 6 ) 

V 2 / n~l 1 + ( nUJ ) 

where a n and b n are the Fourier coefficients of the source h6.1\) . 

Proof. Based on the T-periodicity and the null mean value of g, the proof of 
theorem 15.41 can be straightforwardly extended to prove Y > 0. To prove the second 
inequality, we proceed in three steps. First, as stated in the proof of proposition ^. II 
null forcing gives Y = 0, and hence Y(0) — 0. Secondly, lemma RT2l gives Z(Q) = 0. 
Thirdly, the T-periodic function W = §^ = ^§ satisfies (|4.9I) . The sign of /' and 
/ in (|6.2|) show that the isocline of zero slope (|4.10p is positive or null: 1^ (t) = 
occurs only at points where Z vanishes. The proof of proposition 14.41 therefore still 
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holds here, yielding W(t) > at all t S [0, T]. We therefore obtain W(A) > 0, and 
hence Z(A) > Z(0) — 0, which proves the second inequality. 
Lastly, a Taylor expansion of Y when 4cl gives 



Y(A) = F( )+a|^(0) 



A 2 



A 2 8 2 Y 



2 OA 2 



(0) + O(A 3 ), 



TF(0) + G>(A 3 ). 



Elementary calculations on 



lead to 



W(0) = f (0)Z 2 (0), 

where Z 2 (Q) is given in (|6.5[) . which completes the proof. □ 

In the case of a pure sinusoidal source, a\ = 1 and a n> \ = b n = 0: the estimate 
(|6.6p recovers the estimate (|5.2p . 

7. Conclusion. 

7.1. Physical observables. A set of nondimensionalized parameters was used 
throughout this study. To recover the physical observables, we recall (13.91) and (|6.1[) . 
The results of physical interest are as follows. 

• Theorems 15.41 and 16.31 mean dilatation of the crack 



>0, 



d[u] 

u 13d 



>o, 



T (0) 



P 2 d 



E 



i 



(nn/py 



o 



/3 3 eP 



(7.1) 



where a n and b n are the Fourier coefficients of the normalized source. The 
main aim of the present study was to establish the validity of the inequalities 
in (|7.1|) with any set of parameters. In the case of model 1 (|2.7|) and a 
monochromatic source (|3.8p . the local estimate obtained here matches the 
theoretical expression obtained in our previous study using a perturbation 
analysis method [16) . In the local estimate of (|7.1|) . an improvement of one 
order of accuracy is reached with a purely monochromatic source: see ()5.2j) . 
Lastly, note that the local estimate with model 1 is twice that obtained with 
model 2 flUJ) . 

Theorem 15.61 maximum aperture of the crack (monochromatic source only) 




max 



n 



J TT 



4«o 



2vo 



o 



(3 2 d 



(7.2) 
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• Theorem 15.71 maximum closure of the crack (monochromatic source only) 



-d<d(-rr I <[»]„,„, <o, 

lt , U „ = - 2 f^= + o(-^), ,7.3, 
lim [w] m in = —d. 

It is worth noting that [u] > and [u] m in bounded are purely nonlinear phenomena: 
with the classical linear law (12.4|l . [u] = and [u] m i n is not bounded independently of 

Vq. 

7.2. Acoustic determination of the contact law. One of the applications of 
the present study is the characterization of the crack model, in particular the finite 
compressibility of the crack. This data is crucial in geomechanics and geohydrology, 
where it is linked to the transport of fluids across fractured rocks [3]. If [u], the source, 
and the physical parameters of f^o and fii are known, then the second equation in (|7.ip 
provides a straightforward mean of determining |JF (0)|/d Note that the stiffness K, 
which is involved in j3 via the proposition 13.11 is classically measured using acoustic 
methods |21j . It then suffices to measure [u]. For this purpose, two possible methods 
come to mind. 

The first method consists in measuring the dilatation of the crack mechanically. 
This requires installing two strain gauges around the crack [14] . However, in many 
contexts such as those encountered in geosciences, this is not practicable. The second 
method consists in performing acoustic measurements of the diffracted elastic waves 
far from the crack. The proposition |4T] can be used here, by writing the Fourier series 
of incident, reflected, and transmitted elastic displacements 

ui(x, t) — {cos(f2 1 — ko x) — 1} , 

oo 

u R (x, t) = -Bp + sin n (fit + k x) + R b n cos n (fit + k Q x)} , ^ ^ 

n— 1 

oo 

u T (x, t) = T Q + { T n sinn (fit - fci x) + T b cosn(fii- ki x)} , 

71=1 

where k$ = fl/ci and fci = ft/ci. Definition of [u] implies 

M=T«-i? Q + g. (7.5) 

On the other hand, the causality of the source and the continuity of the stress (|2.3p 
mean that 

rra , r>a PO CO Vo . . 

pi a T +p c R Q = — . (7.6) 



Solving (|73)) and (|7TrJ)) gives 



PO C + Pl Cl px Cl - p C Vq 

H TV ! 

Pici pici Q 

PO CQ + pi Ci Vq 

o "+" z 77- 

p c q 11 
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Purely acoustic measurements of Rq or Tfi therefore give [it] , from which the nonlinear 
crack parameter can be easily determined. 

7.3. Future lines of investigation. The monotonicity of T in (|2.6[) - or equiv- 
alently / in (|3.9[) and (|6.1[) - was the key ingredient used here to prove the properties 
of the crack: the strict increasing of T - or the strict decreasing of / - proves the 
uniqueness and attractivity of the periodic solution. In addition, the strict concav- 
ity of T - or the strict convexity of / - shows that the mean value of the solution 
increases with the forcing parameter. Relaxing these hypotheses may lead to more 
complex situations, and require more sophisticated tools. One thinks for instance to 
the case of hysteretic models of interfaces [TTJ [7] . 

This paper was focused on the normal finite compressibility of a crack. Coupling 
with shear stress is also an important topic to examine [20] . 

Lastly, further studies are also needed on networks of nonlinear cracks [23] . These 
configurations frequently occur in applications. Direct numerical simulations have 
been carried out in [16] , but a rigorous mathematical analysis is still required. 
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